SISSA REF. 193/93/A 
DFPD 93/A/80 
December 1993 



THE THREE-POINT CORRELATION FUNCTION 
OF THE COSMIC MICROWAVE BACKGROUND 
§ : IN INFLATIONARY MODELS 

Q ' Alejandro Gangui^, Francesco Lucchin^, Sabino Matarrese^ and Silvia Mollerach'* 

Xf^ • ^ SISSA - International School for Advanced Studies, 

via Beirut 2 - J^, 34013 Trieste, Italy 
Dipartimento di Astronomia Universitd di Padova, 
^ I vicolo deirOsservatorio 5, 35122 Padova, Italy 

^ Dipartimento di Fisica "Galileo Galilei" Universitd di Padova, 
O ' via Marzolo 8, 35131 Padova, Italy 

^ . ^Theory Division CERN, CH-1211, Geneve 23, Switzerland 

cn ■ 

^! Submitted to The Astrophysical Journal 

O 



Abstract 

We analyze the temperature three-point correlation function and the skewness of the 
Cosmic Microwave Background (CMB), providing general relations in terms of multipole 
coefficients. We then focus on applications to large angular scale anisotropics, such as those 
measured by the COBE DMR, calculating the contribution to these quantities from pri- 
mordial, inflation generated, scalar perturbations, via the Sachs- Wolfe effect. Using the 
techniques of stochastic inflation we are able to provide a universal expression for the en- 
semble averaged three-point function and for the corresponding skewness, which accounts for 
all primordial second-order effects. These general expressions would moreover apply to any 
situation where the bispectrum of the primordial gravitational potential has a hierarchical 
form. Our results are then specialized to a number of relevant models: power-law inflation 
driven by an exponential potential, chaotic inflation with a quartic and quadratic potential 
and a particular case of hybrid inflation. In all these cases non-Gaussian effects are small: 
as an example, the mean skewness is much smaller than the cosmic rms skewness implied by 
a Gaussian temperature fluctuation fleld. 
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1 Introduction 



The recent detection of large-scale Cosmic Microwave Background (CMB) anisotropies by 
COBE (Smoot at al. 1992) has provided us with a unique opportunity to probe the properties 
of the primordial perturbations which affect the CMB temperature distribution through the 
Sachs- Wolfe effect (Sachs & Wolfe 1967). Future years of observation will reduce the level 
of the noise, thus making possible to extract more and more statistical information from 
the data. Besides the classical analysis in terms of the rms temperature fluctuation on the 
smoothing scale of the DMR instrument, the angular two-point correlation function and the 
lowest order multipoles (Smoot et al. 1992; Wright et al. 1992), preliminary results have 
already appeared on a measurement of the three-point function and skewncss on the COBE 
map (Hinshaw et al. 1993; Smoot et al. 1993). The three-point correlation function of CMB 
fluctuations, and its limit at zero angular separation, the skewness, provide an estimate of 
the size of possible deviations from a Gaussian behaviour of the primordial perturbation 
field. 

There are at present two main theories on the origin of primordial perturbations, namely 
that they originated from quantum fluctuations of scalar flelds during an early inflationary 
era or from topological defects produced by a phase transition in the early universe. Present 
observational data have not yet been able to rule out one of them. It was generally thought 
that the main difference between the two models is the statistical distribution of the resulting 
fluctuations: inflationary models were claimed to predict a Gaussian distribution, while 
topological defects a non-Gaussian one. However, it is now clear that, at least for the 
cosmic string case, even if the effect of a single string is clearly non-Gaussian, the combined 
effect of a large number of them results in a quasi- Gaussian distribution. On the other hand, 
it has been realized that inflationary models also predict small deviations from Gaussianity. 
Hence, it is necessary to make more quantitative estimates of the non-Gaussian features 
expected in each model if one intends to use the predicted statistical distribution of CMB 
anisotropies to discriminate between them. 

Besides instrumental noise and contamination by Galactic emission, any analysis of large- 
scale anisotropies is however limited by one more source of indetermination: the so-called 
"cosmic variance" (Abbott & Wise 1984; Smoot et al. 1992), namely the impossibihty of 
making observations in more than one universe, which severely limits our ability to extract 
intrinsic non-Gaussian signals from the data, especially for anisotropy patterns dominated 
by low-order multipoles (Scaramella & Vittorio 1991). This is particularly important when 
considering the temperature three-point correlation function of the COBE map, since most 
present theoretical models of the primordial density fluctuation process predict rather small 
non-Gaussian effects on such large scales, i.e. a quasi-Gaussian anisotropy pattern. Actually, 
the large beam-width of the COBEDMK experiment blurs any information on smaller scales, 
which makes it practically insensitive to the type of non-Gaussian signatures produced by 
topological defects. Cosmic strings predict relevant non-Gaussian features in the CMB at 
angular scales below a few arcminutes (Moessner, Perivolaropoulos &: Brandenberger 1993). 
This makes even more interesting obtaining a quantitative prediction on the non-Gaussian 
CMB features which arc expected to be produced on very large angular scales in inffationary 
models. In fact, every sensible inflation model will produce some small but non-negligible 
non-Gaussian effects, both by the self-interaction of the inflaton fleld, and by the local 
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back-reaction of its self-gravity. It is our interest here to quantify this prediction, with the 
largest possible generality, so that observational results on the CMB angular three-point 
function can be used as a further test on the nature of the primordial perturbation process. 

Unfortunately, we will find that single-field infiation models generally imply mean values 
for the skewness which are well below the cosmic rms skewness of a Gaussian field, which 
confirms and generalizes earlier results based on a simple toy-model (Falk, Rangarajan & 
Srednicki 1993; Srednicki 1993). In this sense, looking at anisotropics on smaller angular 
scales than COBE would probably appear a more promising strategy. Coulson, Pen & 
Turok (1993) and Coulson et al. (1993) consider degree-scale anisotropics produced by non- 
Gaussian field-ordering theories of structure formation. The results of the present paper 
would also apply to these intermediate scales, only provided the window function appropriate 
to the specific experiment is used. 

Falk et al. (1993) first gave a quantitative estimate of the size of non-Gaussian effects 
through a calculation of the three-point CMB correlation function from perturbations gen- 
erated in a simple model, where the infiaton has cubic self-interactions. Luo & Schramm 
(1993) argued that secondary anisotropics could also play a non-negligible role in this re- 
spect. 

Our paper is the first one where the problem is considered in a totally general and 
self-consistent way. We use the stochastic approach to infiation (e.g. Starobinskii 1986; 
Goncharov, Linde & Mukhanov 1987), as a technical tool to self-consistently account for all 
second-order effects in the generation of scalar field fiuctuations during infiation and their 
giving rise to curvature perturbations. We also properly account for the non-linear relation 
between the inflaton fluctuation and the peculiar gravitational potential. Our derivation 
moreover removes a non-realistic restriction to purely de Sitter background made by Falk 
et al. (1993), which is especially important when non-Zel'dovich perturbation spectra are 
considered. 

The plan of the paper is as follows. In Section 2 we give completely general deflnitions 
for the CMB three-point function in a single sky in terms of multipole coefficients. We stress 
that these results actually apply to whichever type of temperature fiuctuations, as well as to 
any angular scale. We also give what we believe is a more operational definition of angular 
three-point function than has been considered so far (e.g. Luo 1993a,b). In the last part 
of this section we anticipate the results obtained from the Sachs-Wolfe effect due to the 
quasi-Gaussian perturbations of the gravitational potential generated during inflation and 
we discuss the effects of the cosmic variance on the possibility of observing a non-zero mean 
value for the skewness. As a by-product of this analysis we compute the "dimensionless" 
rms skewness for a large range of the primordial spectral index of density fiuctuations n. 
All the results reported at this level only depend on three quantities: i) the rms quadrupole 
amplitude Q, which is fixed by observational data; ii) a model-dependent parameter $3 ~ 
a few, related to the amplitude of the skewness; iii) the rotationally invariant multipole 
coefficients Q (suitably normalized to the quadrupole, i = 2), whose specific dependence on 
i is fixed by the spectral index n. Note that infiation is able to produce perturbation spectra 
with essentially all values of n around unity (e.g. MoUerach, Matarrese & Lucchin 1993), 
as it can be needed to match the COBE data (Smoot et al. 1992) with observations on 
smaller scales. The technical derivation in the frame of stochastic inflation of the connected 
two- and three-point function (or its Fourier counterpart, the bispectrum), for the inflaton 
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field first and the local gravitational potential next, is reported in Section 3. Section 4 
deals with the zero-lag limit of the three-point function: the skewness, for which we provide 
a universal inflationary expression, which we then specialize to some of the most popular 
inflaton potentials: exponential (Lucchin & Matarrese (1985)), quartic and quadratic (Linde 
1983, 1985) as well as a simple potential for hybrid inflation (Linde 1993). Section 5 contains 
some general conclusions. 



2 The CMB three— point correlation function 

Our aim here is to compute the CMB temperature three-point correlation function implied 
by possible departures from a Gaussian behaviour of the primordial peculiar gravitational 
potential. As a flrst step we will deflne the connected two- and three-point correlation 
functions of the CMB temperature as measured by a given observer, i.e. on a single microwave 
sky. Let us then deflne the temperature fluctuation fleld ^{x; 7) = {T{x; 7) —To{x))/To{x), 
where x specifles the position of the observer, the unit vector 7 points in a given direction 
from X and To{x) = J ^^T(x; 7) represents the mean temperature of the CMB measured by 
that observer (i.e. the monopole term) . 

The angular two-point correlation function C2(x; a) measured by an observer placed in x 
is the average product of temperature fluctuations in two directions 71 and 72 whose angular 
separation is a; this can be written as 

C2[x-a)= j J ^^d(7i-72-cosa) — (a;;7i) — (a;;72) . (1) 

As well known, in the limit a — > one recovers the CMB variance C2(x) — J ^^[^{x; 7)]^. 
Expanding the temperature fluctuation in spherical harmonics 

A T-l CO £ 

-Tfri^-n) = E E aT{x)WeYr{7), (2) 

e=i m.=-e 

and writing the Dirac delta function as a completeness relation for Legendre polynomials P^, 
we easily arrive at the expression 

C2(f ; «) = ^ E ^^(cos «)Q|(^) >V| , (3) 

where = J2m=-e the previous expressions We represents the window function of 

the speciflc experiment. Setting Wo = Wi = automatically accounts for both monopole 
and dipole subtraction; for £ > 2 one can take We ~ exp 1)(T^ , where a is the 

dispersion of the antenna-beam proflle, which measures the angular response of the detector 
(e.g. Wright et al. 1992). In some cases the quadrupole term is also subtracted from the 
maps (e.g. Smoot et al. 1992); in this case we also set = 0. 

The analogous expression for the angular three-point correlation function is obtained by 
taking the average product of temperature fluctuations in three directions 71, 72 and 73 with 
flxed angular separations a (between 71 and 72), P (between 72 and 73) and 7 (between 71 
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and 73); these angles have to satisfy the obvious inequahties |a — 7I </? <a + 7. One then 
has 

Cs{x; a, P,^) = I —-^ I [ c?cos'i9i25(cos't9i2 — cosa) / c/cos'(923 

J 4:71 Jo 271 J-1 J-l 

/■I AT AT AT 

x5(cosi923 - COS/3) J rfcos^9i35(cos?9i3 - C0S7) — (f;7i) — (x;72) — (x;73), (4) 

where cos = la - Ip and ipi2 is the azimuthal angle of 72 on the plane orthogonal to 71. 
The above relation can be rewritten in a form analogous to Eq.(0), namely 

C,{x- a, /?, 7) = iV(a, Pn)J^I^J ^^(71 " 72 - cos a) 

AT AT AT 

x5(72-73 -cos/?)5(7i-73-cos7) — (f;7i) — (f;72) — (f;73) , (5) 

where N{a,j3,'j) = a/1 — cos^a —cos'^P — cos^7 +2cosacos/?cos7.[| Setting a = /? = 7 = 
in these general expressions one obtains the CMB skewness C^^x) = J -j^[^{x;i)]'^. Also 
useful are the equilateral three-point correlation function (e.g. Falk et al. 1993) and the 
collapsed one (e.g. Hinshaw et al. 1993 and references therein), corresponding to the choices 
a = /? = 7, and a = 7, /? = 0, respectively. Alternative statistical estimators, more suited to 
discriminate bumpy non-Gaussian signatures in noisy data, have been recently introduced 
by Graham et al. (1993). In all the above formulas, full-sky coverage was assumed, for 
simplicity. The effects of partial sky coverage on some of the statistical quantities considered 
here are discussed in detail by Scott, Srednicki & White (1993). 

Following the procedure used above for C2{x; a), we can rewrite the three-point function 
in the form 

C3(f;a,/?,7) = Ar(a,/3,7)f E E a.^a^^a^^^W,, W,,W,3 

^1/2/3 mi, 7712, ma 

xE E ^.(cosa)P.(cos/3)P,(cos7)7^X^'"^7^S^"^^H^™^ (6) 

where the coefficients HT,}^/^""' = J dQ^Y;^'* {^)Y;^' {^)Yf;;^' {^) , which can be easily ex- 
pressed in terms of Clebsch-Gordan coefficients (e.g. Messiah 1976), are only non-zero if 
the indices £i, rrii {i = 1, 2, 3,) fulfill the relations: \ij —ik\ < < \£j+£k\, h + h + = even 
and mi = 1712 + ^3. The collapsed three-point function measured by the observer in x reads 

c,{x-a) = ^ E E PiAcosa)aT,'aT,'aTrm,m,m:,nzt7r (7) 

^1/2/3 mi,m2,m3 

which, for a = 0, gives a useful expression for the skewness. 

So far our expressions have been kept completely general, they would apply to whatever 
source of temperature fluctuations in the sky, through suitable (usually statistical) relations 
for the product of three multipole coefficients appearing in Eqs.(^) and (|^, and to 
whatever angular scale, through the specific choice of window functions W^. However, in 

^ To show that the latter expression is correctly normalized one can expand the delta functions in Legendre 
polynomials and use the relation J^ei"^^ + '^)Pei^)P£iy)Peiz) = f (1 — — — + 2xyz)^^^^. 
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what follows we shall only deal with large angular scale anisotropics originated from primary 
perturbations in the gravitational potential $ on the last scattering surface via the Sachs- 
Wolfe effect (Sachs & Wolfe 1967). In that case ^{x;^) = + ro7), where vq = 2/Ho 
is the horizon distance and Hq the Hubble constant. As noticed by Luo & Schramm (1993), 
secondary anisotropics produced during the mildly non-linear evolution of perturbations 
through the Rees-Sciama effect (Rees & Sciama 1968; Martmez-Gonzalez, Sanz & Silk 1992), 
also imply a non-zero contribution to the connected three-point function, which is still there 
when the underlying primordial linear gravitational potential is Gaussian. Furthermore, 
these authors suggest that this contribution should strongly dominate over the primary 
inflationary one. However, a more detailed analysis (Mollerach et al., in preparation) shows 
that this is not the case on the angular scales probed by the COBEDMR. Therefore, in what 
follows we will only refer to the primary contribution to the three-point function. A recent 
analysis of the three-point function and the skewness of the COBE data has been performed 
by Hinshaw et al. (1993) and Smoot et al. (1993), who find no statistical evidence for non- 
Gaussian signatures beyond those implied by the cosmic variance: while consistent with the 
random-phase hypothesis, these data can only rule out strongly non-Gaussian fluctuations 
on very large scales. On smaller scales, Graham et al. (1993), analyzing the UCSB SP91 
experiment, were able to detect non-Gaussian features, which, however, might also be of 
non-cosmological origin. 

To obtain definite predictions for the statistics described above, one needs to exploit the 
random nature of the multipole coefficients a^. In our case, these coefficients should be 
considered as zero-mean non-Gaussian random variables whose statistics should in principle 
be deduced from that of the gravitational potential. In such a case one should either obtain, 
by Monte-Carlo simulations, different realization of the sky corresponding to different "cos- 
mic observers" , or analytically compute theoretical ensemble expectation values. The latter 
procedure will be followed here to obtain the mean two- and three-point functions. These 
expectation values are of course observer- i.e. x-independent and can only depend upon the 
needed number of angular separations. 

In the frame of the inflationary model, the calculations reported in the following section 
lead to general expressions for the mean two- and three-point functions of the primordial 
gravitational potential, namely 

(<|.(ro7i)$(ro72)) = |^ E(2^ + l)^^(7i " 72)C, (8) 



and 



(<|.(ro7i)$(ro72)$(ro73)) = §^*3 E (2j + 1)(2£+ 1)C,C, 



X [P,-(7i ■ 73)P^(7i ■ 72) + Pj{l2 ■ li)Piil2 ■ 73) + ^i(73 ■ li)Peil3 ■ 72)] , (9) 

where $3 is a model-dependent coefficient. The ^-dependent coefficients are defined by 
(Qj) = (^^+1) Q'^Ci, with Q = {Ql)^^"^ the rnis quadrupole, and are related to the gravitational 
potential power-spectrum P^{k) through Ci = dkk'^P^{k)jf{krQ)/ dkk'^P^{k)il{krQ), 
where is the £-th order spherical Bessel function.^ The rms quadrupole is simply related 

^The possible infrared divergence of this expression for ^ = has no practical effect on observable 
quantities, since the monopole is always removed. 
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to the quantity Qrms-ps defined by Smoot et al. (1992): Q = \/4^Qrms-PslTQ. For the 
scales of interest we can make the approximation Pi^{k) oc /c""^, where n corresponds to the 
primordial index of density fiuctuations (e.g. n = 1 is the Zel'dovich, scale-invariant case), 
in which case (e.g. Bond & Efstathiou 1987; Fabbri, Lucchin & Matarrese 1987) 

Ci = —J — - — — — (10) 
r(£+|-t)r(| + |) 

The equations above allow to compute the angular spectrum, 

(a^'a^^*) = 5ij^i2Smim2-^Cii , (11) 
and the angular bispectrum, 

= if "^sP^.c,, + c,A. + cMnZtTer ■ (12) 

Replacing the latter expression into Eq.(^ we obtain the general form of the mean three- 
point correlation function. Some simplifications occur for the collapsed three-point function, 
for which we obtain 

{C,{a)) = :^3T^<f3 E (2^1 + 1)(2^2 + 1)(24 + l)P,,(cosa)[Q,Q, + CeA, + CM 
xWe,W,,We,J'e,i,e, (13) 

where the coefficients J-'^^t^h = (^tt)^^/ dVL^ J dQ^Pi-^^j ■ ^')P^^{^ ■ ^')P£,^{j ■ 7') may be 
suitably expressed in terms of products of factorials of ii, £2 and £3, using standard relations 

for Clebsch-Gordan coefficients, by noting that J-'kim = (00™) i'^^- Messiah 1976). The 
CMB mean skewness (C^lO)) immediately follows from the above equation for a = 0. A 
plot of the angular dependence of the collapsed three-point function above, normalized to 
the skewness, is reported in Figure 1, for typical values of the spectral index n. 

It was first realized by Scaramella & Vittorio (1991) that detecting a non-zero three-point 
function or skewness for temperature fiuctuations in the sky cannot be directly interpreted as 
a signal for intrinsically non-Gaussian perturbations. In fact, even a Gaussian perturbation 
field has non-zero chance to produce a non-Gaussian sky pattern. This problem is related 
to what is presently known as cosmic variance, and is particularly relevant for fiuctuations 
on large angular scales, i.e. for low-order multipoles of the temperature fiuctuation field. 
One way to quantify this effect is through the rms skewness of a Gaussian field {C^{0))causs- 
It is easy to find 

{ClmGauss = 3 f_^d COS a{C2{a)f . (14) 

Scaramella & Vittorio (1991) and, more recently, Srednicki (1993), focused on the most 
popular case of a scale-invariant spectrum, n = 1 and the COBE DMR window function, 
corresponding to a Gaussian with dispersion a = 3°2 (e.g. Wright et al. 1992). We will 
be interested here in the same quantity, but for various values of n. In Figure 2 we have 
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plotted the normalized rms skewness (C*! (0))^''f„j,^/(C2(0))^''^ as a function of the spectral 
index, both including and removing the quadrupole: in both cases this ratio is in the range 
0.1 — 0.3 for interesting values of n. The values obtained for n = 1, both with and without 
the i = 2 contribution, are identical to those given by Srednicki (1993), who adopted the 
same smoothing angle (the slightly different definition of window function cannot affect our 
dimensionless skewness ratios). As a rough criterion, we can conclude that, in order to 
detect primordial non-Gaussian signatures, {C^{0)) must be at least of the same order as 
{C^{0))Gauss- other hand, as we will see below, single-scalar-field inflationary models 

generally lead to skewness ratios (C3(0))/(C2(0))^''^ < 10~^, so that their non-Gaussian 
features cannot be distinguished from the cosmic rms skewness. In this sense the quasi- 
Gaussian inflationary predictions for the CMB anisotropies are in full agreement with the 
recent analysis of the three-point function and the skewness from COBE data (Hinshaw et 
al. 1993; Smoot et al. 1993). 



3 Stochastic inflation and the statistics of the gravita- 
tional potential 

Now we will show the validity of Eq.(^ for the three-point correlation function of the 
gravitational potential due to perturbations produced during an inflationary epoch in the 
early universe. This calculation will provide the primordial power-spectrum, reflecting into 
the £ dependence of the Q coefficients, as well as a general expression for the factors Q 
and $3 above. In order to take into account all the effects contributing to a non-vanishing 
primordial three-point correlation function of $, we will perform the computation in two 
steps. In Section 3.1 we compute the three-point function for the inflaton field perturbation 
50. The most convenient way to perform this calculation is in the frame of the stochastic 
approach to inflation (Starobinskii 1986, Goncharov et al. 1987), which naturally takes 
into account all the multiplicative effects in the inflaton dynamics that are responsible for 
the non-Gaussian features. Then, in Section 3.2 we compute the extra-contribution to the 
three-point function of the gravitational potential that arises due to the non-linear relation 
between $ and 5(f). This effect has been previously noticed by Barrow & Coles (1990) and 
Yi & Vishniac (1993). 



3.1 The inflaton bispectrum 

To study the dynamics of the inflaton, we will apply the stochastic approach. This is based 
on defining a coarse-grained inflaton field 0(x, a), obtained by suitable smoothing of the 
original quantum field over a scale larger than the Hubble radius size, whose dynamics is 
described by a multiplicative Langevin-type equation. This is obtained by adding to the 
classical equation of motion a Gaussian noise term whose amplitude is fixed by the rms 
fluctuation of the scalar field at Hubble radius crossing, 

= -x:^ + m,^s a) (15) 
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where V{4>) is the inflaton potential, primes denote differentiation with respect to and 
K = vovrG = \^87r /uip, with mp the Planck mass. The Hubble parameter here should be 
consistently calculated from the local energy density of the coarse-grained inflaton. The 
noise term rj has zero mean and autocorrelation function (e.g. MoUerach et al. 1991) 



{ri{x, a)ri{x, a')) = hjo{qs{a)\x — x'\)6{a — a' 



(16) 



The use of the time variable a = In {a/a^) in this equation has been motivated by Starobin- 
skii (1986), who noticed that a accounts for the possible time dependence of the Hubble 
parameter. This is particularly relevant when general, i.e. non-de Sitter, inflation is stud- 
ied. We also defined the coarse-grained domain size through the comoving wave-number 
qs{a) = eH{a)a{a), with e a number smaller than one, H{a) = H{(j)ci{a)), with (pd{o:) the 
homogeneous classical solution of the Langevin equation (i.e., that obtained with the noise 
term "switched" off). Finally the scale-factor a{a) is obtained by integration of H{a). 

The stochastic dynamics of the coarse-grained field within a single coarse-graining do- 
main (i.e., for X = x') can be studied in terms of the Fokker-Planck equation for the proba- 
bility distribution function of (f). In our case, instead, since we are interested also in spatial 
correlations of the field, we will solve directly the Langevin equation above. To the aim of 
computing the three-point function of 0, a second-order perturbative expansion around the 
classical solution is enough. We will require that the potential ^((/>) is a smooth function of 
its argument, which translates into requiring well defined values for the steepness of the po- 
tential X(a) = X(0c«(a)) = 'mpV {(pel) IV {(pel) and its derivatives (Turner 1993) throughout 
the range of relevant scales. Apart from this requirement we keep the analysis general; only 
at the end we will apply our results to some specific inflationary potentials. We first expand 
V^(0) around (pcu up to second order in (50(x, a) = 0(x, a) — 0d(«), (i-e. up to order ^), 
\^(0) = V{(t)ci) + V'{(t)ci)^4' + |^"(0c/)'^0^ + ■ ■ ■. Replacing this into the Langevin equation 
we obtain 



d6(f){x, a) 
da 



A{a)5(j){x, a) + B{a)5(l)'^{x, a) + [Di{a) + D2{a)S(j){x, a)] r]{x, a) , (17) 



where we have used d(f)ci/da = —V'{(f)ci)/i^'^V{(j)ci)- In Eq.(|T^ we defined 



A = -^X'-B 
Stt 



167r 



27r 



Let us now split the field perturbation as 6(j) = 6(t>i + of 0{h^^'^) and 0{h), respectively. 
We can also define a rescaled variable fj = {H{a) /2Ts:)r]. We then find 



i)i(af, a) = X(a) / da'X {a')ri{x,a') 
Jo 



(19) 



MX, a 



X(a) 



da' 



Bia' 



X(a 



—501 {x, a') H S(j)i{x, a')r]{x, a') 

2m- 



Let us now calculate the connected three-point correlation function of 
non-vanishing contribution {0{h'^)) reads 



(20) 



The lowest order 



{S(j){xi, ai)6(f>{x2, a2)S(j){x3, as)) 
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as 

da' 



Jo 



^^'^j-{6(f)i{xi, ai)6(f)i{x3, a')) {6(j)i{x2, a2)S(t>i{x3, a')) + [xi^ X2] 



c?a'[(50i(fi,ai)50i(x3,a'))('^0i(^2,a2)r?(f3,«')) + [^i^^2]] 

Imp Jo 

+ 2x2 terms . (21) 
The term proportional to X{a3)/2mp in the r.h.s. of this equation can be recast in the form 

^ ^"1 da'X-\a'y-^{6MSi,c^i)SMS3,a'))joiqsia')\x2-X3\) (22) 



2mp Jo (27r) 

where we defined amin = min[a3,a;2]. We need now to compute the 6(j) auto-correlation 
function. To this aim, recalling that a is the time when the perturbation wavelength a/q 
equals the size of the coarse-graining domain, we change the integration variable in Eq.(p!9D) 



from a to g = ^9>nV{a)/^ q^e'^ / H^mp. The subscript * denotes quantities evaluated at the 
time when we start to solve the Langevin equation; this is chosen in such a way that the 
patch of the universe where we live is homogeneous on a scale slightly above our present 
horizon (see e.g. the discussion by MoUerach et al. 1991). We then find 

1 Mrmniquq') XiqAXiq') 

(50i(xi, ai)5(j)i{x3, " )) = ^ dqq ^(^)^(^jY(^^o(gl^i - ^sl) (23) 

where we defined P{q) = \q^^H^{q) where a = a{q). Using this we can rewrite Eq. ([2^) as 
X{q^)X{q2)X{q3) f d^q' f d\ P{q')P{q) 



2mp J (27r)3y {2tt)^ X^{q) 

xW{q'; g™„(g3, q2))W{q; g^,„(gi, q')y'ii^i--s)^i<fi^2-.-s) (34) 

where we defined the filter function W{q; qi) = 6(g — g*) — 6(g — g^) (9 is the Heaviside 
function). 

A similar analysis can be performed for the term proportional to B{a) in Eq.(^). We 

get 

X{q,)X{q2)X{q3)l^ -jB{q)X{q)J J x^i^q") X-(q"') 

xW(g"'; g™„(g2, q'))W{q"- q^^Qu q'))d'^'<^^-^^)e''^"-^^^-^^^ . (25) 

So far we have been working in configuration space. In order to obtain the gravitational po- 
tential at Hubble radius crossing during the Friedmann era, it is convenient to Fourier trans- 
form the coarse-grained infiaton fluctuation. One has 6(j){x, <y{q)) = {27t)^^ J d^k 6(f){k)<d{q — 
k)e^^'^, where 6(f){k) denotes the Fourier transform of the full scalar field at the time a{q). 
This follows from the fact that at the time a{q) the only modes k that contribute to the 
coarse-grained variable are those which have already left the inflationary horizon, namely 
k < q. We can then obtain the Fourier transform 6(f){k,a{q)) = 6(j){k)Q{q — k), and, in the 
limit k q^ (or equivalently q —>■ k"^, that is, when we consider the horizon-crossing time 
of the given scale), we simply have 6(f){k, a{k)) = 6(j){k). In other words, at the time a{q) the 
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Fourier transform of the coarse-grained variable coincides with that of the full field. Using 
these results in Eq.(0) we finally obtain the inflaton bispectrum through 



((50(A;i, a{h))6<P{k2, a{k2))6(j){k,, aik^))) = {2ny6'{k, + k^ + k,)P{k2)P{k, 



X 



X(A;i 



X\k2)Q{k2 - ks) + X^{k,)eiks - k^ 



X{k2)X{k,) 

+{ki^k2} + {h- 



2m I 



+ 2 



da' 



(26) 



3.2 The three— point function of the gravitational potential 

We want now to compute the three-point function of the peculiar gravitational potential, 
($(ro7i)$(ro72)$('"o73))- The approximate constancy of the gauge-invariant quantity C out- 
side the horizon (e.g. Bardeen, Steinhardt & Turner 1983) allows to obtain the gravitational 
potential during the matter-dominated era, given the value of 5(j) during inflation. During 
inflation one has C(x, a) ~ — 50(x, a) / {dip/ da), which is usually interpreted as a linear rela- 
tion between C, and 50. However, when calculating the three-point function of $ one cannot 
disregard the second-order effects coming from the fluctuations of d(t)/da. 
Recalling that 

(90 mpV'((j)) mp , . , 

one gets 

This equation is expressed in configuration space; the Fourier transform of the first term 
inside the brackets is just 50(fc, a). For the second term we get a convolution of the type 
/ d'^p6(f){'p)6(f){k — p). For the scales of interest one has $(fc) ~ —3C{k)/5. Then, adding the 
two above contributions and evaluating the expression at the horizon crossing time we get 



$(A;): 



S^Aa{k))-^^^^l-0^J<P{p,a{k))5<P{k^^^ 



5mpX{a{k)) 
From this equation we calculate 

mhMk^Mh)) = x{k^)xl^^^^^ 

' ^ ^ ^ ' ■{6(f){p)6(t){ki-p)6(j){k2)6(f){ks)) 



(29) 



X^{k^)X{k2)X{ks) J (27r)3 
+{h^k2} + {h^h} (30) 

Using the fact that, at horizon crossing, (50(A;i)50(fc2)) = {2TT)^P{ki)5^{ki + k2), we can 
write 

mh)^k2)) = {2Txf5\k, + k2)f\a{k,))P{k,) , (31) 
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with f{a{k)) = 2A'K/5mpX{a{k)). Finally, using Eqs. 



H) and ( 



i^) we find 



X 



x\k2)eik, - k,) + x\k^)Q{k^ - k. 



2m I 



X\k2)X\k^) 

2X\k^)X{k^)X{k^^ 



X\k^) 



dq' 



+2 / -^B{q')X{q') \ + {k^^k^} + {h ^ k^} 



q 



(32) 



We are interested in considering perturbation modes that left the horizon about 60 e-foldings 
before the end of the inflationary epoch. In the explicit examples below X{k) turns out to 
be a slowly varying function of k. We therefore approximate X{k) ~ Xqq in what follows. 
Then, to the lowest non-vanishing order, the three-point correlation function is 



{^{k{)^{k2mk:,)) 



1 

/eo 



X, 



60 



2m I 



2^60 



+ 



60 



X, 



60 



'^eo dq 



B{a{q))X{a{q)) 



27r) V(A;i + k2 + fcg) [P^{kr) P^^k^) + P^{k2)P^{k^) + P^{k^) P^{kr)] 



(33) 



with /eo = 247r/5mpX6o. Using the expression for P{q) we obtain the power-spectrum for the 
pecuhar gravitational potential P$(/i;) = 2'k'^ f^^k''^ H"^ {a{k:)) / Atx'^ ^ \flQHlQk~^{k/k^Y^^ , 
up to possible presence of small logarithmic corrections. The primordial spectral index n 
is related to the inflationary parameters by the approximate relation n ~ 1 — (XIq/Stt) + 
[mpX'^Ql Ait) (Turner 1993). Now we can obtain the two- and three-point correlation func- 
tions for $ in configuration space by inverse Fourier transforming Eq.(PTD and Eg. (P3[) . We 
then get Eq.(||) and Eq.(P|) with 





r(3- 


n)T{l 


+ f) 


bmlXio 


r (2- 


t)l 


2 

r 


'9 n\ 
•-2 2' 



(34) 



and 



1 

/60 



60 



2X' 



60 



2m £ 



X 



+ 



60 



^=60 dq 
Xgo Jk, q 



B{a{q))X{a{q)) 



(35) 



In order to obtain Eq.(|^) and Eq.(^, we had to use the definition of The one given 
in Section 2 differs from the results obtained here because of the window function W{k), 
appearing in the inflationary expressions; this can be however neglected if we account for 
the oscillating behaviour of ji for large arguments and for the fact that ji — for small 
arguments (that helps in cancelling the lower divergence). 



4 The CMB skewness 

In what follows we will restrict our analysis to the mean skewness as given by Eq.(|T3p, for 
a = 0. In the numerical calculation below we will assume that the higher multipoles are 
weighted by a 7°5 FWHM beam, resulting in cr = 3°2 (e.g. Wright et al. 1992). To normalize 
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our predictions we will consider the rms quadrupole obtained by Seljak & Bertscliinger (1993) 
through a Maximum-Likelihood analysis, namely Qrms-ps = (15.7±2.6) exp [0.46(1 — 'ri)]yU-ft' 
(see also Scaramella & Vittorio 1993 and Smoot et al. 1993). Q For the models with n < 1, 
considered below, this expression should be multiplied by an extra factor [(3 — n)/{lA — 
12n)]^/^ to account for the effective decrease in the estimated value of Qrms-PS due to the 
contribution of gravitational waves (Lucchin, Matarrese & MoUerach 1992). For the mean 
temperature we take the FIRAS determination Tq = 2.726 ± O.OIK (Mather et al. 1993). 

To estimate the amplitude of the non-Gaussian character of the fluctuations we will 
consider the "dimensionless" skewness Si = {Cs{0)) / {C2{0))^^'^ ■ Alternatively, if we want 
our results to be independent of the normalization, we may also define the ratio ^2 = 
((73(0))/ (C2(0))^, as suggested by the hierarchical aspect of our expression for the skewness. 
We obtain 



60 



^60 



2^3/2 (^) 



and 



^2 



15 
167r 



X 



60 



Ami 



XL 



60 



-^60 



+ g 



To(n) 



(36) 



(37) 



with Q = AmpX^Q J^^°{dq/q)B{a{q))X{a{q)) and where we also defined the spectral index- 
dependent geometrical factor 

^^^^^ ^ lEi,,e,,e,mi + l){2e2+l){2i3+l)[C,A.^^^^ ^3g^ 



where the exponent p in the denominator takes values 3/2 and 2 for Si and S2, respectively. 
The numerical factors Ip{n) in Eq.(|38D are plotted in Figure 3 for different values of the 
primordial index. 



4.1 Inflationary models 

Let us now specialize our general expressions to some simple inflationary models. 

Exponential potential 

Let us first consider power-law inflation driven by the exponential potential V{(j)) = 
Voexp (— Ak0), with A < (Lucchin & Matarrese 1985). In this case the power-spectrum 
is an exact power-law with n = 1 — 2A^/(2 — A^). We note in passing that the right spec- 
tral dependence of the perturbations can be recovered using the above stochastic approach 

■^Note that this value of Qrms-PS assumes a muhivariate Gaussian distribution function, accounting 
for both the signal and the noise; while in principle one should repeat the analysis consistently with the 
assumed statistics of the temperature perturbations, the quasi-Gaussian nature of our fluctuation field 
allows to extrapolate this Maximum-Likelihood estimate without sensible corrections. 
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(Mollerach et al. 1991). For this model we find X 
A = B = 0. We then have 



— vSttA, whose constant value implies 



Si 



3A if 60 

4 mp 



3/2 



n] 



r(3-n) 


r(| + f) 


r (2 - f ) 


;'r(|-t) 



1/2 



So 



15 



(39) 



The COBE results constrain the amplitude of i^eo- For the case n = 0.8 we have Hgo/mp = 
1.8 X 10~^. This gives Si = 9.7 x 10~^ and iSi = 1.1 x 10"^, without and with the quadrupole 
contribution respectively, while 1S2 = 1.3 in both cases. 

Quartic potential 

Consider now the potential ^X(f)^ for chaotic inflation (Linde 1983, 1985). For this model 
X = 4mp/0 and therefore mpB ~ rripX" -C X. Thus we can take ^ ~ in the coefficients 
of Eqs. (|37|) and (|36D for the dimensionless skewness parameters. By using the slow-roll 

where the 
ends when 



^60 



1 + Hk/ho) 



solution for the infiaton we may write X ~ Amp 

logarithmic correction to the scale invariant power-spectrum is small. Inflation 
the infiaton takes the value (pend - 
spectral index is n ~ 1. We find 



2/377 mp implying 



^60 




30 / mp" 



60/7r mp. In this case the 



(40) 



COBE constrains the value of A to be A ~ 1.4 x 10 We get Si = 5.2 x 10 ^ and 
Si = 6.0 X 10^^ without and with the quadrupole contribution respectively, while ^2 = 0.5. 

Quadratic potential 

Another simple potential for chaotic inflation is |m^(/)^ (Linde 1983, 1985). In this case 
X ~ 2mp06o^ [1 + (2/K^0go) ln(A;//c6o)], 06o — y^SO/V mp and n ^ 1. We find 

3 m^ 



Si 



2y^ mp 



3/2 1 



So 




Ml) 



(41) 



with m^/mp ~ 1.1 x 10 ^ from the COBE normalization. We get Si = 3.9 x 10 ^ and 
Si = 4.5 X 10~^ without and with the quadrupole contribution respectively, while ^2 = 0.4. 

Hybrid inflation model 

Finally, let us consider a model of hybrid inflation recently proposed by Linde (1993). 
Inflation happens in this model during the slow-roll evolution of the inflaton in the effective 
potential V^(0) = Vq + \m^(l?, where Vq = M^/4 is a cosmological-constant-like term. At a 
given time, when the inflaton fleld takes the value = M, its coupling with a second scalar 
fleld triggers a second-order phase transition of the latter (whose vacuum energy density is 
responsible for the cosmological constant term), which makes inflation end. An interesting 
prediction of this model is that the spectral index of density fluctuations for wavelengths 
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which left the horizon while Vq > m?(f)^/2 is larger than unity, namely n ~ 1 + 2m? / hi^Vo 
(Mollerach et al. 1993). For this potential we find X = mpiv?(f)/VQ, implying Q = Q. We 
then obtain 



Si 



73 060 

2 mpM2 




r(3- 


n)V[l 


+ f) 


r (2- 


t)l 


2 

r 


'9 n\ 
-2 2' 



1/2 



(42) 



So 



15 iv^piri^cf)] 



60 



TC 



MS 



60. 



(43) 



The slow-roll solution for the infiaton gives 06o = Mexp[60 m^/ft^Vo]. Choosing a value 
n = 1.1, the COBE results constrain the free parameters to be M ^ 1.3 x 10~'^mp and 
m ~ IQ'^nip. In this case we find Si = —1.1 x 10^^ and Si = —1.3 x 10^^ without and 
with the quadrupole contribution, while S2 = —1.6. Note that these results suggest some 
correlation between the sign of the skewness and the value of the spectral index n. 



5 Conclusions 

The results reported in the previous section seem to preclude any chance of actually obtaining 
observable non-Gaussian signals at least in the frame of infiation, unless one resorts to more 
complicated multiple-field models (e.g. Allen, Grinstein & Wise 1987; Kofman et al. 1991; 
Salopek 1992). However, it should be stressed that all the results reported in Section 2 would 
apply to any large-scale anisotropy where the temperature fiuctuation can be obtained by a 
local perturbative calculation (this is not the case, for instance, for secondary anisotropics, 
which are to be ascribed to an integrated effect): under these conditions the hierarchical 
form of Eq.(^, where the bispectrum is a sum of products of two power-spectra, holds. 

One may also wonder whether the general analysis of the CMB skewness generated in the 
frame of the inflationary model may provide further constraints on the model parameters. We 
did not address this issue in the present paper, where we only considered some simple realistic 
models. In these relevant constraint can be obtained. Both for the exponential 

potential and for the polynomial ones the skewness can be estimated by taking just the flrst 
term in the brackets of Eqs.(^) and (0). An upper limit on its magnitude is then obtained 
by requiring that they give rise to an accelerated universe expansion (i.e. inflation), which 
provides a constraint on the steepness of the potential, X"^ < 247r. This corresponds to 
Si < 10-" and ^2 < 20. 

Notwithstanding the non-zero value of the skewness, any actual detection of this signal is 
hardly distinguishable from the cosmic variance noise in which it is embedded: the intrinsic 
limitation induced by our impossibility of making measurements in more than one universe. 
In fact, the overall coefficient of Si is generically much smaller than the dimensionless rms 
skewness calculated from an underlying Gaussian density fleld. Removing the quadrupole 
contribution to reduce the cosmic rms skewness, as recently suggested by Luo (1993b), does 
not change this main conclusion, because the predicted mean skewness would also be reduced 
by a comparable factor. Things get even worse if we take into account that no sampling 
of the CMB sky is complete, not even the COBE one, where a cut for galactic latitudes 
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\b\ < 20° is required. As shown by Scott et al. (1993), this "sample variance" contributes 
in hiding the signaL Therefore, any possible skewness detection on the COBE scale is most 
probably due to these statistical effects than to any primordial non-Gaussian feature. 
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Figure captions 



Figure 1: Collapsed three-point function (with and without the quadrupole contribution), 
normalized to the skewness, as a function of the angular separation a (left panel). 
The different curves refer to three typical values of the primordial spectral index: 
n — 0.8, 1, 1.2. The symbols in the left panel are the same as in the right one, which 
shows a vertical expansion of the same plot. 

Figure 2: Normahzed rms skewness of a Gaussian temperature fluctuation field as a func- 
tion of the spectral index n, both including and removing the quadrupole contribution. 

Figure 3: Geometrical factors Xp(n) as a function of n. Crosses are for J3/2(n), where the 
quadrupole term is removed. Diamonds correspond to 23/2(^1) including the quadrupole 
contribution. Open boxes correspond to X2(n). For the latter case the difference in 
including or removing the quadrupole contribution is O(10~^), which is not perceivable 
in this graph. 
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